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Abstract 

Hyperbolic monopole motion is studied for well separated monopoles. It is shown that the motion 
of a hyperbolic monopole in the presence of one or more fixed monopoles is equivalent to geodesic 
motion on a particular submanifold of the full moduli space. The metric on this submanifold is found 
to be a generalisation of the multi-centre Taub-NUT metric introduced by LeBrun. The one centre 
case is analysed in detail as a special case of a class of systems admitting a conserved Runge-Lenz 
vector. The two centre problem is also considered. An integrable classical string motion is exhibited. 
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1 Introduction 



One of the crowning achievements of classical mechanics was the demonstration by Newton that Kepler's 
laws of planetary motion followed from an inverse square law for gravity. The publication of Newton's 
Philosophiae Naturalis Principia Mathematica is considered by some to mark the beginning of modern 
physics. Mathematically, as well as historically, this problem is of particular interest. Bertrand's theorem 
iTfl shows that the only central forces that result in closed orbits for all bound particles are the inverse- 
square law and Hooke's law. The reason for these closed orbits in both cases is the existence of 'hidden' 
or 'dynamical' symmetries. These are symmetries which only appear when one considers the full phase 
space of the problem. In both cases the dynamical symmetry group is larger than the naive 50(3) arising 
from the geometrical rotational symmetry of the central force problem. For the inverse-square force, 
on whose generalisations we shall focus, the dynamical symmetry group is 50(4) in the case of bound 
states and 50(3, 1) in the case of scattering states. The group is enlarged due to the existence of an 
extra conserved quantity, known as the Runge-Lenz vector although it pre-dates both Runge and Lenz. 
Certainly Laplace and Hamilton both knew of this conserved quantity but the identity of the first person 
to construct it is unclear. It is claimed in [2 1 that the credit for this belongs to Jakob Hermann and Johann 
I. Bernoulli in 1710, 89 years before Laplace and 135 years before Hamilton. Following common usage 
however, we shall refer to the 'Runge-Lenz' vector. Systems which admit such a conserved vector are of 
great interest, and have been the subject of much research over the years. 

The 'hidden' symmetries also appear in the quantum mechanical problem. The system may be quan- 
tised by the usual process of replacing Poisson brackets with commutators and in this way the energy 
levels and degeneracies of the hydrogen atom may be found. The degeneracy is greater than expected 
since the states of a given energy level fit into an irreducible representation of 50(4) rather than of 
50(3). This explains why the energy depends only on the principle quantum number and not on the total 
angular momentum as one might expect. 

The conservation of the Runge-Lenz vector also implies that the orbits of a particle moving in a New- 
tonian potential are conic sections. Thus the mechanics of Newton is elegantly related to the geometry 
of Euclid. In the early 19th century Lobachevsky and Bolyai independently discovered that other ge- 
ometries are possible which violate Euclid's parallel postulate: the spherical and hyperbolic geometries. 
An interesting question is to what extent the special properties of the Kepler problem in flat space carry 
over to these new geometries. The generalization of the Kepler problem on spaces of constant curvature 
was studied by Lipschitz Q and Killing and later rediscovered and extended to the quantum case 
by Schrodinger Q and Higgs J6) 3 . It was found that a Runge-Lenz vector may be constructed which is 
again conserved. Unlike in Euclidean space, the dynamical symmetry algebra is not a Lie algebra. The 
problem on hyperbolic 3-space H 3 is closely related to that on 5 3 and most results are valid for both 
when expressed in terms of the curvature. The problem in hyperbolic space is simpler in a sense, since 
global topological restrictions prevent, for example, a single point charge from existing on the sphere. 

The proof of Bertrand assumes that the force on the particle depends only on its distance from the 
origin. A particle moving in a magnetic field will experience velocity dependent forces and so violates 
the assumptions of the theorem. The simplest magnetic field that one may consider is that due to a 
magnetic monopole at the origin. It was discovered by Zwanziger [8] and independently by Mcintosh 
and Cisneros @ that a charged particle moving in the field of a monopole at the origin has closed 

3 For a review of the history of this problem see Q 
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orbits if the potential takes the form of a Coulomb potential together with an inverse square term. This 
problem, known as the MIC-Kepler or MICZ-Kepler problem can also be generalised to the sphere and 
the hyperbolic space ifTTl ITOl [T2l . There is once again a hidden symmetry responsible for the closed 
orbits, which is a generalisation of the Runge-Lenz vector. For a review of these systems, see lfP3l . 

A final system which exhibits a conserved Runge-Lenz vector is that of two well separated BPS 
monopoles in E 3 . These are solitonic solutions to 5 , t/(2)-Yang-Mills which behave like particles at 
large separations. For slow motions, it is possible to describe the time evolution of the field in terms 
of a geodesic motion on the space of static configurations, the moduli space lTT4l . In the case of two 
monopoles scattering, the full moduli space metric was determined by Atiyah and Hitchin 1151 . In the 
limit when the monopoles are well separated, this metric approaches the Taub-NUT metric with a negative 
mass parameter. The problem of geodesic motion in Taub-NUT is in many ways similar to the motion of 
a particle in a Newtonian potential. In particular, there is once again a conserved vector of Runge-Lenz 
type which extends the symmetries beyond the manifest geometric symmetries lfT6l [P71 . The orbits as 
a consequence are conic sections and the energy levels of the associated quantum system have a greater 
degeneracy than might be expected. 

The question which we consider in this paper is to what extent the results concerning well separated 
monopoles in flat space can be extended to hyperbolic space. In section El we show that for SU (2) 
monopoles in H 3 one can make progress by assuming that some of the monopoles are fixed at given 
positions. The motion is then described by a geodesic motion on a space whose metric is a generalisation 
of the Taub-NUT metric where one considers a circle bundle over H 3 rather than E 3 . Such metrics 
were first introduced by LeBrun lfl~8l . In section|3]we discuss some of the geometric properties of these 
metrics. In section |U we consider the classical mechanics of a quite general Lagrangian admitting a 
Runge-Lenz vector which includes as special cases all of the systems mentioned above, in particular the 
LeBrun metrics. Also included in this section is a discussion of a particular classical string motion in the 
LeBrun metrics which admits a Runge-Lenz vector. In section|5]we consider the quantum mechanics and 
derive the energy levels and scattering amplitude in the context of the general Lagrangian. In section|^lwe 
consider the problem with two fixed centres and show that this is integrable both in classical and quantum 
mechanics, generalising a known result about the integrability of the ionised hydrogen molecule. Finally 
we collect for reference in the appendix some useful information about hyperbolic space. 

2 Hyperbolic Monopoles 

We are interested in SU (2) monopoles on hyperbolic space. These are defined as follows ll"H?l. Let H 3 be 
hyperbolic 3-space of constant curvature —1, with metric h. For a principle SU(2) bundle P — > H 3 , let 
(A, $) be a pair consisting of a connection A on P and a Higgs field $ £ T(H 3 , g P ), where gp denotes 
the associated bundle of Lie-algebras P x 5(7(2) su{2). The pair (A, $) is a magnetic monopole of mass 
M £ M>o and charge k £ Z> if it satisfies the Bogomol'nyi equation 



d A ® = - *„ F A 



(2.1) 



and the Prasad-Sommerfeld boundary conditions 
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Here S 2 is a sphere of geodesic radius r about a fixed point in H 3 and F A is the curvature of the connection 
A. Hyperbolic monopoles were introduced by Atiyah ll20ll . who constructed them from S" 1 -invariant 
instantons on S 4 . We define the moduli space to be the space of solutions of (12.11 12.21) modulo gauge 
transformations, i.e. the space of physically distinct solutions. Atiyah also showed that the moduli space 
of hyperbolic monopoles of charge k can be naturally identified with the space of rational functions of 
the form 

,/ v ai^ 1 + d2Z k - 2 + ... + a k 

f( z > = — m k-i i n — (2 - 3) 

z k + o\z* 1 + . . . + b k 

where the numerator and the denominator have no common factor. Thus the moduli space is a manifold 
of dimension (Ak — 1), provided k > 1. In fact this was only shown for the case when 2M is an integer. 
Following [19] we shall assume that the moduli space is a manifold of this dimension for all values of M. 
The moduli space is often taken to be enlarged by a circle factor. This can be defined by fixing a direction 
in H 3 , say X\, choosing the gauge so that Ax — and then only considering gauge transformations which 
tend to the identity as the hyperbolic distance along x\ tends to infinity. From now on we shall consider 
this enlarged moduli space, which has dimension Ak. 

For BPS monopoles in flat space, it is possible to find a metric on the moduli space for well separated 
monopoles by treating them as point particles carrying scalar, electric and magnetic charges ll2T1 . We 
denote by M. k the moduli space of fc-monopoles. The A;-fold covering of M. k , M. k splits as a metric 
product E3- 

M k = M° k x S 1 x E 3 . (2.4) 

The E 3 component physically represents the centre of motion of the system, which separates from the 
relative motions governed by the metric on Ai^. The S 1 factor is an overall phase, which may be thought 
of as a total electric charge which is conserved. The fact that the centre of mass motion can be split off 
in this fashion is due to the Galilean invariance of the theory, which appears as a low velocity limit of 
the full Poincare invariance of the S , f/(2)-Yang-Mills-Higgs theory. In the case of two well separated 
monopoles Mt® is Euclidean Taub-NUT with a negative mass parameter. 

For monopoles in H 3 , the dimension of the moduli space is known to be the same as that for flat 
space. We call the moduli space of fc-monopoles on hyperbolic space Mk and its fc-fold covering M k . 
Little is currently known about the geometry of this moduli space, although it is known that the L 2 metric 
one defines in the flat case is infinite. What this means for monopole scattering is unclear. It is certainly 
plausible that the slow motions are still described by a moduli space geodesic motion. We give support 
to this idea below by showing that this is the case in the unphysical situation where k — 1 monopoles are 
at prescribed locations. We shall assume that the metric on the space of such configurations arises as the 
restriction of a metric describing the interaction of k free monopoles, however our results do not rest on 
this being the case. 

We may separate off an S* 1 factor, corresponding to the phase at infinity, so we might expect a metric 
on M k to have a metric product decomposition analogous to (12.41 ) of the form 

M k = A? xS'x H 3 , (2.5) 

however there is no analogue of the Galilei group for hyperbolic space. Translational invariance guar- 
antees a conserved total momentum, but since there are no boost symmetries, this cannot be split off as 
a separate centre of mass motion. In general if one considers two particles in H 3 interacting by a force 
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which depends only on their separation, one finds that the relative motion depends on the total momen- 
tum of the system ll23l . ll24ll . Thus, it is not unreasonable to expect that the moduli space metric will not 
split up in this fashion. This means that if we wish to study the scattering of two hyperbolic monopoles, 
we will be studying geodesic motion with respect to a 7-dimensional moduli space metric. 

In order to make this problem more tractable, we shall consider a simplified situation where one or 
more monopoles are fixed at given locations. This is unphysical in the sense of SU (2) monopoles since 
for well separated monopoles, the mass of each is not a free parameter, it is determined by the other 
charges and so we cannot consider an infinitely heavy monopole in order to fix it in one place. Although 
this situation does not represent a genuine motion of hyperbolic monopoles, it is nevertheless of interest 
as it gives an insight into the interaction of two hyperbolic monopoles. We find that the equations of 
motion may be put into the form of geodesic motion on a generalisation of the multi-centre metrics, first 
written down by LeBrun |[T%ll . 

A similar situation is considered by Nash ll25ll . who has investigated singular monopoles on H 3 . Here 
one allows the solution of the Bogomol'nyi equations to be singular at some points but the nature of 
the singularity is fixed. Nash shows using twistor theoretic methods that the natural metric for a charge 
1 singular monopole is indeed a LeBrun metric. Unlike the metrics we find below which generalize the 
singular negative mass Taub-NUT metrics, Nash finds metrics which generalize the everywhere regular 
positive mass Taub-NUT metrics. 

The LeBrun metrics share a lot of the properties of the multi-centre metrics, including integrable 
geodesic equations for the one and two centre cases. This allows progress to be made analytically and 
indeed makes the metrics worthy of study in their own right. It is also worth noting that the multi-centre 
metrics may be found as hyper-Kahler quotients of the asymptotic k monopole moduli space metrics 
(26|. We may conjecture that these LeBrun metrics may arise from the full hyperbolic monopole moduli 
space metric by some analogous construction. 

2.1 Motion of a test monopole 

In order to simplify the problem of hyperbolic monopole scattering, let us consider n monopoles fixed in 
H 3 at points Pj, I — 1 . . . n. Following Manton ETI . we treat these monopoles as point particles carrying 
an electric, magnetic and scalar charge. This is possible since outside the core, the 577(2) fields abelianise 
as in flat space. If the radius of the hyperbolic space is large compared to the core of each monopole, 
the mass and scalar charge will be given in terms of the other charges by the Euclidean expressions. We 
wish to consider a test particle moving in the field of n identical monopoles, so we take the electric and 
magnetic charges to be q and g 4 respectively for all of the monopoles. The scalar charge of a monopole 
with magnetic charge g and electric charge q is (g 2 + q 2 ) 2 . We wish to interpret the electric charge as a 
momentum, and since in the moduli space approximation all momenta must be small, we assume that q 
is small. 

As in the case of flat space, the electric and magnetic fields E, B are one-form fields on H 3 , defined 
in terms of the field strength 2-form by: 

F — dt A E + * h B. (2.6) 

4 Unlike in the previous section, we do not take units such that the magnetic charge is an integer. 
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The field-strength 2-form, away from Pj, satisfies the vacuum Maxwell equations: 



dF = 0, d*F = 0. (2.7) 
In the above equations, * is the Hodge duality operator for the metric 

ds 2 = -dt 2 + h, (2.8) 

where h is the metric on H 3 , with Hodge duality operator -k h . Clearly, from (12.71) there is a symmetry 
under F — > *F. This is the usual duality, well known in flat space, that for a vacuum solution one can 
replace (E, B) with (B, —E) and Maxwell's equations continue to hold. Because of this duality it is 
possible to locally find two gauge potentials, A and A, which are related to F by 

F = dA, F — *dA. (2.9) 

We will require both these potentials, since a monopole moving in this electromagnetic field will couple 
to both of them. The potential and dual potential at a point P due to a monopole at rest at a point Pj is 
given by: 

A = -^-Vdt+^-uj 

4.7T 4.7T 

A = -—Vdt-—uj 1 (2.10) 
4tt 4% 

where V is a function on H 3 and uj is a one-form on H 3 , which are defined by: 

V = coth D(P, Pj) - 1, * h dw = dV, (2.11) 

with D(P, Pi) the hyperbolic distance between P and Pj. Since V satisfies Laplace's equation on H 3 , 
d *h dV = 0, uj may be determined up to a choice of gauge (which we will return to later). Using these 
equations, it is possible to check that dA = *dA, which implies Maxwell's equations (12.71) for the 2- 
form F. Since Maxwell's equations are linear, it is straightforward to write down the potential and dual 
potential at a point P due to all n of the monopoles at points Pj. This is still given by (12.101 ). but with V 
and uj defined by 



V = coth D ( p ' p i) 



n 
i=i 

dV = * h duj. (2.12) 

There is also a scalar field in this theory, the Higgs field, $. The value of $ at the point P due to the 
monopoles at Pj is given by 



(g 2 + q 



2- 



2 



$ = A V, (2.13) 

4tt 

with V as in equation (12.121) . We now consider another monopole at the point P, with velocity vector 
v E TpH 3 . This monopole has magnetic charge g, electric charge q' and mass m. The Lagrangian for 
the motion of this monopole in a field with potential A and dual potential A and with Higgs field $ is: 

L = [-m+ (q 2 + q' 2 )^}(-g(v,v))^ + q'{A,v) + g{A,v), (2.14) 
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Where g is the metric given by (I2.8I ). v = (1, v) is the 4-velocity of the particle and (•, •) is the usual 
contraction between vectors and one-forms. This is the form of the Lagrangian argued by Manton in 
||2T1 . generalised to hyperbolic space. 

As mentioned above, we are interested in slow motions of the monopoles, such as might be described 
by a moduli space metric approach. Thus, we will assume that v, q and q' are all small and expand to 
quadratic order in these quantities. We define q = q' — q and take the coordinates of P to be r\ i = 1 , 2, 3, 
so that v l = r\ We also add a constant term m — mq 2 /2g 2 to the Lagrangian, which will not affect the 
equations of motion. The new Lagrangian after the expansion is: 

-m V r hijr J + — ujir -m V q . (2.15) 

2 8tt J 3 An g 2 \2 8n J H 

We note that if we can interpret q as a conserved momentum, this Lagrangian will be quadratic in time 
derivatives and so may be interpreted as a geodesic Lagrangian. In order to do this, let us consider a 
circle bundle, E, over H 3 \ {Pi}, which is the configuration space of our problem. We assume a metric 
form: 

ds 2 = Ah + B(dr + Wf (2.16) 

with h the metric on H 3 and W a connection on the bundle. This gives a geodesic Lagrangian of the 
form: 

L = ^Ar'h.jH + \b{+ + rW,) 2 . (2.17) 
Since r is a cyclic coordinate, there is a conserved quantity which we will identify as a multiple of q: 

q = K B(f + r i W i ). (2.18) 

As (12.171) is a Lagrangian, we cannot simply use (12.181) to eliminate f, but must instead use Routh's 
procedure (see for example ll27l0 which gives a new effective Lagrangian with the same equations of 
motion as (I2.17I) . This new Lagrangian is given by: 



1 a ■ la 2 

-Af'haf 3 + -r l Wi 1 

2 k 2k 2 B 

Comparing this with (12.151) . we see that they are the same if we identify 



U = -.\rh, :i H + IrWj - -4yv (2-19) 



A = m-fv, Wt ^ = ~ 2 A. (2.20) 

47r 47r B g l 

We fix the value of k to remove the Misner string singularities due to u, and this gives k = A-k/ g. In 
order to clean up the form of the Lagrangian, we take units where m = An and g = Airp and remove 
an overall factor of An. We thus have that the slow motion of the monopole is equivalent to a geodesic 
motion on a space with metric at P given by: 

ds 2 = Vh + V- 1 (dr + uf (2.21) 

Where h is the metric on H 3 and we redefine V so that 

y=l+p^n-^cothD(P,Pj)V (2.22) 
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Finally, we still have that uj is defined in terms of V by the relation 

* h duj = dV. (2.23) 

This only defines uj up to an exact form on H 3 , however we can see from ( 12.211) that a change of gauge 
for uj simply corresponds to a coordinate transformation on E. 

Thus we find that the motion of a monopole in H 3 in the presence of n fixed monopoles may be 
described by geodesic motion on a hyperbolic multi-centre metric given by (12.221) . (12.23b . Geometrically 
we may think of this as the asymptotic metric on a submanifold within the full moduli space metric. This 
however will not be a geodesic submanifold since there is no mechanism to fix n of the monopoles while 
allowing one to move. 

3 The LeBrun metrics 

The metric ( 12.211 ) was considered by LeBrun with p = —1/2 in |[T%ll . For this choice of p, the metric 
is everywhere regular, the apparent singularities at Pj are nuts, i.e. smooth fixed points of the U(l) 
action generated by the Killing vector J^. This can be seen by considering an expansion in a small 
neighbourhood of one of the Pj where the curvature of the hyperbolic space can be neglected. The 
metric then looks like that of Taub-NUT. In the case we consider above with positive p, the metric will 
become singular near the fixed monopoles where V becomes negative. Since we wish to consider an 
approximation where all the monopoles are well separated, we can ignore this singularity and we expect 
that it is smoothed out in the full moduli space metric. This is precisely what happens in the case of 
monopoles in flat space, where the full Atiyah-Hitchin metric is regular, while the negative mass Taub- 
NUT is not. 

LeBrun was interested in finding a metric on CP 2 # • • • #CP 2 which was Kahler and scalar flat. He 
showed that in the conformal class of metrics of which (12.211) is a representative, there is a whole 2- 
sphere's worth of different metrics which are scalar flat and Kahler, one corresponding to each point at 
infinity of H 3 . More explicitly, he showed that if V and u are as in (12.221) . (12.231) and q is any horospherical 
height function, then the metric 

g = q^Vh + V^idr + uj) 2 ) (3.1) 

is Kahler with scalar curvature zero. Here a horospherical height function is a function on H 3 whose 
restriction to some geodesic is the exponential of the affine parameter and which is constant on the 
forward directed horospheres orthogonal to this geodesic. The horospheres are discussed in section dA.21 ). 
Since this construction picks out a point on the boundary of H 3 arbitrarily, there is a two-sphere's worth 
of conformal factors which make the metric (12.211) Kahler. These metrics have self-dual Weyl tensor and 
so are sometimes referred to as half-conformally-flat. 

LeBrun in addition showed that the metric g represents a zero- scalar-curvature, axisymmetric, asymp- 
totically flat Kahler metric on C 2 which has been blown up at n points situated along a straight complex 
line, i.e. n points have been replaced with CPi. This aspect of the metrics was considered in ll28l where 
C 2 with an arbitrary number of points blown up was considered as a spacetime foam for conformal 
supergravity. 

Since the LeBrun metrics are all conformally scalar flat, an interesting question is whether there is 
a conformal factor such that the metric is locally Einstein. This was studied by Pedersen and Tod ll2"91 
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and they found that a metric of the form ( 12.211) is locally conformally Einstein if and only if (V, to) is a 
spherically symmetric monopole. In other words only the one-centre metrics are conformally Einstein. 
A Kahler form K is both closed and covariantly constant 

dK = 0, V/i = 0. (3.2) 

Thus it necessarily satisfies Yano's equation: 

dK = 3VK (3.3) 

Returning to the conformal representative given by (12.211) then, we find that this metric admits a 2- 
sphere's worth of conformal Yano tensors which satisfy the conformal Yano equation: 

K XK , a + K aK , x = (g aX K v K - >v + g^K^) . (3.4) 

The existence of such conformal Yano tensors permits the construction of symmetry operators for the 
massless Dirac equation on this space QUI . 

We show below that in the case of the one- and two-centre metrics there are additional higher rank 
symmetries which echo those of the multi-Taub-NUT metrics. The one-centre metric has Killing vectors 
which generate the SU (2) x U(l) isometry group which can be made manifest by writing the metric in 
a Bianchi-IX standard form. In addition, there are 3 rank two Killing tensors which transform as a vector 
under the SU(2) action and as a singlet under the U(l) action. These we interpret as a generalisation of 
the Runge-Lenz vector. The two-centre metric has the reduced isometry group of U(l) x U(l) and in 
addition admits a rank 2 Killing tensor. This generalises the extra conserved quantity associated with the 
motion of a particle in the gravitational field of two fixed centres. 

We see then that the LeBrun metrics have many attractive geometric properties which are closely 
related to those of the multi-Taub-NUT metrics and they are worthy of closer study. 



4 Classical mechanics of the one-centre problem 
4.1 Geodesies on Self-Dual Taub-NUT 

We begin by recalling some results about the geodesies of Self-Dual Taub-NUT lfT6ll3Tll32ll . Taub-NUT 
is a four dimensional Hyper- Kahler manifold, with metric: 



ds 2 



a 



3- 



(4.1) 



where <jj are the usual left invariant forms on SU (2). The geodesic Lagrangian for negative mass Taub- 
NUT (we take p = — 2, since changing p only changes the overall scale of L) is 



r • r + 4 1 



(ip + cos 



(4.2) 
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There is a conserved quantity q associated to the ignorable coordinate ip and also a conserved energy E. 
These are given by 

? = 4(l-H) (^ + 008^), E =\{ 1 - 2 ^) ( i-2+ (|) 2 )- (43) 

The conserved angular momenta are given by: 

J = (l - -J r x r + qr , (4.4) 

where f is a unit vector in the r direction. In addition to these conserved quantities, there is a extra 
conserved vector, 

K = ( 1 - - j f x J+ \2E- -q 2 ) f. (4.5) 



rj V 2 . 

These conserved quantities allow all of the geodesies to be found. 

As one would expect, the conserved quantities are related to geometric symmetries of Taub-NUT. 
Taub-NUT has biaxial Bianchi-IX symmetry, so the group of isometries is 577(2) x U(l). The SU(2) 
factor gives rise to the conserved vector J which transforms as a triplet under the SU (2) action. The 
U(l) factor gives rise to the conserved SU(2) singlet q. The origins of the vector K are less obvious. 
Since K is quadratic in the momenta, it arises from a triplet of rank 2 Killing tensors on Taub-NUT, each 
satisfying Killing's equation: 

= 0- (4-6) 

These Killing tensors are best understood in terms of the Kahler structures of Taub-NUT. Since Taub- 
NUT is Hyper- Kahler, there are three independent complex structures F^ v . In addition, Taub-NUT admits 
a further rank 2 Yano tensor, i.e. an anti-symmetric tensor obeying Yano's equation: 

Y^o) = 0. (4.7) 

The Kahler structures, being covariantly constant, also satisfy Yano's equation trivially. It was shown in 
P2"l that the Killing tensors arise as the product of the complex structures with the extra Yano tensor: 



4.2 A class of systems admitting a Runge-Lenz vector 

We wish now to consider the generalisation of (14.11 ) which describes the motion of a hyperbolic monopole 
about a fixed monopole, as we derived in section |2~T1 This corresponds to geodesic motion on a manifold 
with metric (I2.21I) . We will find that many of the integrability properties of geodesic motion in Taub- 
NUT persist. Rather than start with the Lagrangian for geodesic motion, we instead consider a more 
general Lagrangian and impose conservation of a Runge-Lenz vector. We find that within the class of 
Lagrangians of this form admitting a Runge-Lenz vector is the Lagrangian describing geodesic motion 
on (I2.21I ). We shall consider the Lagrangian: 



V(\t\) 



r • r 



(1 



I |2N 2 



V(|r|) 



W{\r\ 



(4.9) 
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We note that the kinetic part is given by a metric of the form (I2.21I ). with h the metric on H 3 in Beltrami 
coordinates (see section lA. 1.3ft . We write this in spherical coordinates, and assume a form for uj such 
that the Lagrangian (14.91 ) admits a rotational SU(2) symmetry. Then: 



2 



V(r) 



j. 2 



iy 2 j 2 ^2 



) 2 + sin 2 i 



T + p COS 



Since r is cyclic, the conjugate momentum p T = e is conserved. This is given by: 

e = V~ (f + pcosi 

There is also a conserved energy, E: 



E 



V(r) 



? 2 + sin 2 ^ 2 )+e 2 



+ W(r) 



(4.10) 



(4.11) 



(4.12) 



1 — r 2 (1 — r" 

If we define a new variable ^ by r = then the Lagrangian (14.101) is manifestly SU (2) invariant. Thus 
there is a conserved angular momentum vector, which we can write as: 

V(r) 



-r x r 



1 — r* 
r x N + epr , 



epr 



where we have defined for convenience a new vector 

N = V(r) 



1 



(4.13) 



(4.14) 



We would like to find the condition that the Lagrangian ( 14.91 ) admits a Runge-Lenz vector in addition 
to the constants of the motion above. By analogy with the Taub-NUT case, we posit a further conserved 
quantity of the form: 

K = NxJ + #r. (4.15) 

Using the Hamiltonian formalism it is a matter of straightforward, if tedious, calculation to check that K 
is a constant of the motion if and only if the functions V(r), W(r) and the constant g satisfy the equation 



-e 2 V(r) 2 + (W(r) — E)V(r) 



fj 



e 2 p 2 



C, 



(4.16) 



2 ' v/ 7 ' v/ r 2r 2 

where C is an arbitrary constant. Now V(r) and W(r) cannot depend on E, as they are given functions 
in the Lagrangian (14.9ft . However g and C are free to depend on E (see for example (14.5ft ). To find the 
most general functions satisfying (14.161) for some g(E), C(E) we differentiate twice with respect to E 
and find: 



V{r) 



9\E) 
r 

g"(E) 



C\E), 
-C"{E). 



From (14.171 ) we can write 

g = - a E + (3, c=-kE + j. 
Substituting these back into (14.171) and (14.161) we prove the following proposition: 



(4.17) 
(4.18) 
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Proposition 4.1. The most general functions W(r), V(r) such that the Lagrangian \4.9\) admits a Runge- 
Lenz vector of the form (14.751 ) are given by: 

V(r)=k + i, ^(f + g+^+H)"'-^^). (4.19) 

From now on, we assume that V(r) and W(r) take this form. 

4.2.1 Some Special Cases 

We can find some special cases of the Lagrangian (14.91 ) by choosing particular values for the constants: 

• If we take k = 1, a = and e = 0, we have the Lagrangian for the Kepler problem on hyperbolic 
space, H 3 . 

• If we take k = 1, a = 0, but allow e to be non-zero, we get the Lagrangian for the MICZ-Kepler 
problem in hyperbolic space lfT2ll . It will be useful for later to note that (formally at least) we can 
also obtain this system by setting a = p and letting p — > while requiring that ep = fi remains 
constant. 

• If we take a = =pp, (3 = —kae 2 , 7 = e 2 k 2 /2 and k — 1 ± p, we have the Lagrangian for pure 
geodesic motion on the manifold with metric (I2.21I) . 

In specifying the metric on H 3 , we have everywhere assumed that the space has radius 1. Rather 
than carry an extra parameter through the calculations, it is easier to replace the hyperbolic radius R 
by dimensional analysis in any equations. Doing this, we can consider allowing R to tend to infinity. 
The metric on H 3 approaches that on E 3 , so we might expect to recover known results for flat space. 
In particular, with the special parameter choices above, we may expect to recover results for the Kepler 
problem, MICZ-Kepler problem and the problem of geodesic motion on Taub-NUT. 

4.2.2 Poisson Algebra of Conserved Quantities 

We have 8 conserved quantities, E, e, J, K, which fall into two singlet and two triplet representations of 
SU (2) respectively. We can calculate the non-vanishing Poisson brackets and after some algebra we find 

{K^Kj} = -2e ijk J k (Ek + J 2 - e 2 p 2 - 7) . (4.20) 

Thus the Poisson algebra of the conserved quantities is not a Lie algebra, since the bracket of two com- 
ponents of K is a cubic polynomial in the conserved quantities. Algebras of this kind are known as finite 
jy-algebras. They have been extensively studied in J33j. This particular algebra is some deformation of 
the so(4) or so(3, 1) algebra, depending on the sign of Ek — e 2 p 2 — 7, which leaves the so(3) subalgebra 
undeformed. These algebras are considered in P41I351 . 
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4.2.3 The shape of the orbits 

The existence of extra constants of the motion for the system defined by (14.91) and (14.191) means that it is 
possible to determine the shapes of the orbits entirely from the conserved quantities. In the Beltrami ball 
model (see lA.1.3l we find that the orbits are given by conic sections. 
From the definition of J, (14.131) . we can immediately see that 

J ■ r = epr. (4.21) 

This is true for any functions V(r) and W(r) and it means that the motion lies on a cone centred at the 
origin, with axis parallel to J and half angle 9 satisfying cos6> = ep/ |J|. In the case where ep = 0, we 
find that the motion is in a plane orthogonal to J. The SU(2) symmetry therefore allows us to reduce the 
three dimensional problem to motion on a two dimensional surface. In general however, the orbits will 
not be closed and the motion may be chaotic. 

When we have a conserved Runge-Lenz vector however, any bound orbits must be closed. To see 
this, we consider K • r. 

K-r = r-NxJ + gr, 
= J • r x N + gr, 

= J (J - epr) + — J • r. (4.22) 

ep 

Where we have used the definitions (14.131) . (14.151) and the equation (14 .21b . Collecting these terms, we 
find that for ep ^ 

9 t \ t2 „2J2 



K - — J = J - eV- (4-23) 

ep J 

This is the equation of a plane orthogonal to K — — J which in general does not pass through the origin. 
Thus the orbits lie on the intersection of a cone and a plane and are hence conic sections. 

If we take ep = 0, then r ■ J = and K ■ J = 0. The motion is then in a plane perpendicular to J 
and we may take the coordinates on this plane to be r and 9, the angle that r makes with the vector K. 
Taking r ■ K then yields the equation 

IKI \ J 2 

1 - ^cosfl = . (4.24) 

9)9 

This is the general form for a conic section in the plane. Thus, we have shown that for the Lagrangian 
defined by (14.91 ) and ( 14.191 ) the general orbit is a conic section in Beltrami coordinates. 

It should be noted that for ep ^ 0, the plane in which the motion lies does not include the origin. 
Let us consider an incoming particle, scattered by some infinitely heavy particle at the origin, such that 
the motion is described by our Lagrangian. The two body system defines a plane which includes the 
locations of both the particles and the velocity vector of the moving particle. In the Kepler problem in 
flat space, the motion of the particles is fixed in this plane, however here the plane of motion receives a 
twist as the incoming particle is scattered. This phenomenon is noted in [16] and is typical of monopole 
scattering. 

It is possible to consider all of the special cases mentioned above and we find agreement with the 
results in the literature, as reviewed in the introduction. We can also include the hyperbolic radius R , 
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and in the limit R — > oo we recover the flat space results for the Kepler and MICZ-Kepler problems, 
and also the results for geodesic motion on self-dual Taub-NUT. 

In conclusion, we have found that for the general Lagrangian admitting a Rung-Lenz vector given 
above, the orbits are always conic sections in an appropriate set of coordinates. 



4.2.4 The Hodograph 



In a paper of 1 847 [|36| Hamilton introduced the hodograph, defined to be the curve traced out by the 
velocity vector, thought of as a position vector from some fixed origin, as a particle moves along its orbit. 
He showed that the Newtonian law of attraction gave rise to a hodograph which is a circle, displaced 
from the origin, and conversely that this was the only central force law to give a circular hodograph. 

We shall define the hodograph in H 3 to be the curve swept out by the vector N, thought of as a 
position vector originating at O. By differentiating the condition (14.231) with respect to time, and using 
the definition of N, (14.141) . it can be seen that the hodograph lies in the plane 

N • (gj - epK) = 0. (4.25) 

We may construct coordinates such that N lies in the x\ — £2 plane, then by squaring the definition of K 
it is possible to show that the hodograph is an ellipse, with centre at 

J 2 K 2 - e 2 p 2 g 

x K (4.26) 



J 2 |J x K| 



and eccentricity given by 



/ J\g 2 -K 2 ) \ / e 2 p 2 f g 2 -K 2 Y 2 2 g 2 - K 2 
" =eP { 9 - 9 (J 2 -e 2 P 2 ) ) l^b^vj - 2eP l 2 ~^ 2 V 2+K 1 • ^ 

We see that in the case where e or p vanishes, the hodograph is circular. In the case where a also vanishes, 
we find that the Kepler problem does indeed give a circular hodograph. For a/Owe find a more general 
Lagrangian with a circular hodograph, however this corresponds to modifying the kinetic term, so there 
is no contradiction with Hamilton's result in flat space that the unique central force law with a circular 
hodograph is the Newtonian interaction. 



4.3 The Hamilton- Jacobi Method 

An extremely powerful method for finding solutions to a problem in Hamiltonian mechanics is the 
Hamilton- Jacobi method. For a Hamiltonian H(q\pj), we seek solutions S(q i ) to the Hamilton-Jacobi 
equation: 

dS 

HW, W ) = E. (4.28) 

If a general solution to the Hamilton-Jacobi equation can be found then it is possible to solve the system 
completely, i.e. give the positions and momenta at a time t, as a function of t and the initial positions and 
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momenta at some time t . On the other hand, a particular solution of (14.281 ) defines a coherent family of 
orbits whose momenta at each point are given by: 

» = f (**> 

Both these approaches to solutions of the Hamilton-Jacobi equation are discussed by Synge 071 . We will 
show that the Hamilton-Jacobi equation separates, and then consider a particular solution corresponding 
to Rutherford scattering in hyperbolic space. 

The Hamiltonian for the Lagrangian ( 14.91) . after replacing p T with e, may be written 

H = \ (V 1 ( Pi - eui) h ij ( Pj - ecuj) + Ve 2 ) + W. (4.30) 

Since we have already solved the problem of motion governed by this Lagrangian in Beltrami coor- 
dinates, we shall instead consider the solution in pseudoparabolic coordinates (see section lA. 1.51) . This 
will allow us to construct surfaces of constant action for particles scattered by the origin. This was done 
in the case of a Coulomb potential in flat space by Rowe Il3"&ll3"^1. 

We will require some results from section |A. 1 .5[ Taking cylindrical polars p, z, <p on the Beltrami 
ball, we can define the pseudoparabolic coordinates fi, v to be given by: 



/i 2 — v 2 \xv 



(4.31) 



2 + /i 2 -z/ 2 ' 2 + /i 2 -z/ 2 ' 

Where 

0</i<oo, < v < 1. (4.32) 
In these coordinates the metric on H 3 takes the form: 

dfj? dv 2 



ds 2 = -, » 2 + u2 



v 2 )(l+H 2 ) 



+ ~ o + 



1 — V 2 jJL 2 + V 2 



(4.33) 



We again assume that uj = p cos 9d(p. From the equations (14.311 ) and the standard relations between 
spherical and cylindrical polar coordinates, it is a straightforward matter to express r and cos# in terms 
of [i and v. We then have all of the necessary information to form the Hamilton-Jacobi equation (14.28b . 
We find that the equation separates if W and V are given by (I4.19I ). with the additional condition that 
a = p, i.e. we require that equation ( 12.231) holds. Explicitly, we can write: 

S = M(fi) + N(v) + m0 + er, (4.34) 

then M(/i) and N(v) satisfy the ordinary differential equations: 

2 ,(dM\ 2 (m + ep) 2 2E{k-p) + 2(3-2 1 -e 2 p 2 



" ' - -K. (4.35) 

K is an arbitrary constant which arises from the separation of variables. It gives an extra constant of the 
motion related to the Runge-Lenz vector. In the case of pure geodesic motion, we find that K is quadratic 



dN\ 2 (m - ep) 2 2E(k + p) + 2/5 + 2 7 + e 2 p 2 
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in the momenta, and so corresponds to a second rank Killing tensor of the metric (12.21b . We have picked 
out a particular direction in constructing our coordinates, so it is clear that there is a rank 2 Killing tensor 
associated to each point on the sphere at infinity. 

The differential equations (14.351) can be solved analytically, but for illustration we shall consider a 
simpler case. As was previously noted, if we take the limit p — > 0, with e = and k = 1, we recover 
the Kepler (or Coulomb) problem on hyperbolic space. We may set 7 = without loss of generality, 
since this simply corresponds to shifting the zero point of the energy. We do not seek to find the general 
solution of (14.28b . Following Rowe, we would like to find a solution which describes a family of orbits, 
initially parallel to the positive z-axis 5 and which are scattered by the Coulomb centre. 

Since the particles move initially parallel to the z-axis, they have no angular momentum about the 
2:- axis, so we set m = 0. We define 

k = y / 2(E + 0) and A = y/2(E - 0). (4.36) 

The equations for M and N then simplify considerably, and we find 

(1+ ">Ur) =k ~tt?' =—>- K - (437) 

We shall consider the case of an attractive Coulomb centre, i.e. [3 < 0. 6 We take K = A 2 in order to 
simplify the second equation. This is necessary in order to have the correct v dependence as fi — > 00, i.e. 
far from the origin along the positive z-axis. Integrating then gives: 

N = ^log(l-„ 2 ) 



M = ^ log V ? + " ' + KIX ) + A log ( ±A„ + VV(1 + ^)-^ (4.38) 




There are in principle two choices of sign to be made, one for k and one for A, however it can be seen 
that k — > — k does not change the equations (14.381) . If we consider the solution with the upper signs, as 
P -> we find that 

S(», log(l + /i 2 )(l - v 2 ) + /(k, A) (4.39) 

The function / is singular in the limit k — > A, however S is only defined up to an additive constant, so 
it can be removed. The level sets of the function (1 + yU 2 )(l — u 2 ) are known as horospheres and play a 
special role in scattering in hyperbolic space. They are the hyperbolic equivalent of the plane wave-fronts 
in flat space, see section lA~2l 

The solution (14.391) represents a family of free particles which move parallel to the z-axis. If we allow 
j3 to be non-zero, (14.381) represents a family of particles initially moving parallel to the z-axis but which 
are scattered by a Coulomb centre at the origin. We need to consider the choice of sign however. In fact, 
we need both branches of (14.381) to describe the scattering process. Since H 3 \ {0} is not topologically 

5 We say that two geodesies are parallel if they meet on <9H 3 . When we say a geodesic is parallel to the z-axis, we implicitly 
mean that it meets the positive z-axis on the boundary of H 3 . Two geodesies which do not meet, even on the boundary of H 3 
are ultra-parallel. 

6 The repulsive case follows in a similar fashion to the attractive case 
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Figure 1: A plot showing the surfaces of con- 
stant action with E = 10, (3 = — 1 for the whole 

ofH 3 



Figure 2: A plot showing the surfaces of con- 
stant action with E = 10, (3 = — lin the region 
around the origin 



trivial, we may allow multiple valued solutions to (I4.28I) . Alternatively, since both branches are equal on 
fi = 0, we can follow Rowe and define S as a single valued function on two copies of hyperbolic space, 
communicating through the negative z-axis. A particle which starts initially parallel to the positive z-axis 
and is bent through the negative z-axis will pass from one copy of H 3 to the other. 

Taking /i — > oo, we see that the Coulomb action approaches the free action, but with some distortion 
of the plane wave. This is typical of the Coulomb potential, even in flat space. 

In order to see physically what these solutions mean, it is convenient to plot the surfaces of constant 
action. For the Hamiltonian above, with the parameters as defined, Hamilton's equations give 

r j = h lj Pj (4.40) 

thus the orbits are everywhere normal to the surfaces of constant action. If we plot these surfaces in the 
Poincare model of H 3 (see lA.1.21 . then the orbits will cross the surfaces of constant S at right angles in 
the Euclidean sense. Figures [T] and |2] show plots in Poincare coordinates of the curves S = const, in a 
plane containing the z-axis. Both branches of S have been plotted on the same graph, and we see that 
they match along the branch cut. We also see that one branch represents incoming waves and the other 
gives scattered wavefronts. 

4.4 A classical string motion 

In 114*01 classical Nambu strings were studied on a multi-monopole background of the form 

g AB dx A dx B = -dt 2 + V^idx 4 + Uidx 1 ) 2 + Vh ij dx l dx i (4.41) 



17 



with A, B = 0, . . . , 4, i,j = 1, 2, 3. Here hij = 5ij is the metric on E 3 , x 4 is a compact direction such 
that < x 4 < 2itRk, and V and oj take the usual form: 

k 

V = 1 + — V -, — - — r , duo = * h dV. (4.42) 

2 / — ' x — xJ 

i=i 

It was found that a classical string winding the x 4 direction m times behaved like a relativistic particle 
moving geodesically in an effective 4-metric. Furthermore, this motion may also be described by a 
classical motion where the strings are attracted to the monopoles by a Newtonian inverse square force. 
Thus the relativistic string motion in the one or two monopole case has all of the symmetry of the classical 
Kepler and two gravitating centre problems. 

Classical superstrings wound on the x 4 direction of (14.411) were considered in iHTI in the context of 
Type II string theory compactified on T 6 . It was shown that since the surfaces of constant r correspond 
to squashed 3-spheres which are simply connected, it was possible to unwind a string whilst keeping it 
arbitrarily far from a central monopole. Since the winding number of a string is a conserved charge, the 
H-electric charge, it was argued that the monopole should be able to carry this charge and the zero mode 
which carries it was identified. 

An interesting question is whether there is an analogous behaviour for winding strings in the LeBrun 
metrics. Since the LeBrun metrics are not Ricci flat, it is not straightforward to exhibit a superstring 
solution whose string metric is LeBrun. We shall therefore consider the case of classical Nambu strings 
moving in a LeBrun background. We shall see that the string motion may be thought of as particle motion 
in H 3 , with each monopole attracting the string with a force given by the hyperbolic generalisation of the 
Newtonian attraction. 

We assume that we have a classical string moving in a background with metric (14.411k where now 
is the metric on H 3 and V takes instead the form 

k 

V = 1 + — J2(cothD(P,Pi)-l), du = * h dV. (4.43) 

i=i 

The equations of motion are derived from the Nambu-Goto action, which may be written 



S = —, I ( det -g A B d * Ad f h ) du l du 2 . {AAA) 



2ira' J \ du a du h 

The motion of the string is specified by giving x A as a function of the world-sheet coordinates u a = (a, r) . 
We shall consider a closed string, so that the embedding functions must satisfy the periodicity conditions 

x A ( a , t) = x A (a + 2tt, r). (4.45) 

The simplest way to satisfy these equations is to consider a string which winds m times around the 
compact 'internal' direction, x 4 . The ansatz we shall make is that the functions x A (a, r) take the form: 

x 4 = maR K , x a = x a (r), (4.46) 

where a = 0, . . . , 3. By finding the equations of motion for x A (a, t) and imposing the conditions (14.461) 
it is possible to show that the equations of motion for x a (r) take the form of geodesic equations for a 
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relativistic particle of mass ji = mR k /a' in an effective metric: 



g a pdx a dx 13 = —V 1 dt 2 + hijdx l dxK 



(4.47) 



Thus we may study the string motion by studying the Lagrangian 



L — hj a x *x? i 



d_ 



(4.48) 



We take r to be an affine parameter so that L = — e/i 2 , with e = —1, 0, +1 corresponding to spacelike, 
null and timelike geodesies respectively, t is clearly cyclic, so may be eliminated in favour of its conjugate 
momentum E and we arrive at the energy conservation equation: 



Thus the string coordinates x a behave like a particle of mass m moving in HI 3 in a potential (e — V)E 2 /2. 
From (14.431) it is clear that this potential is that due to k attractive particles, whose force on the string is 
given by the hyperbolic Coulomb interaction. 

We have studied the case of one or two fixed Coulomb centres in H 3 as special cases of the systems 
considered elsewhere. In particular, for the one centre case, there is a conserved angular momentum and 
also a conserved Runge-Lenz vector. In Beltrami coordinates, the path of the string will be a conic sec- 
tion. Thus, if the string initially has some angular momentum about the monopole, they can never collide. 
If the string however starts with no angular momentum, it will eventually reach the point where the x 4 di- 
mension, and hence the string, has zero proper length. It seems plausible then that these monopoles catal- 
yse the annihilation of these winding strings, as was conjectured in the case of Kaluza-Klein monopoles 
over flat space. 

5 Quantum Mechanics of the one-centre problem 

It was argued in llT6l that in a low energy regime the quantum behaviour of two BPS monopoles may be 
approximated by considering a wavefunction which obeys the Schrodinger equation on the two-monopole 
moduli space, whose metric is the Atiyah-Hitchin metric. Since in this case the motion is geodesic, the 
Hamiltonian is proportional to the covariant Laplacian constructed from the Atiyah-Hitchin metric. If we 
are interested in quantum states where the amplitude for the particles to be close together is small, we 
can find approximate results by considering the self-dual Taub-NUT metric as an approximation to the 
Atiyah-Hitchin metric for well separated monopoles. 

We are interested in the behaviour of well separated monopoles in hyperbolic space. As was noted in 
sectional the absence of boost symmetries means that the case of two particles interacting via a Coulomb 
interaction does not reduce to a centre of mass motion together with a fixed centre Coulomb problem 
ll24l . Thus, we expect that the problem of two monopoles in hyperbolic space will require the solution 
of Schrodinger's equation on a 7-dimensional manifold. In order to make analytic progress, we instead 
consider the problem of the motion of a test monopole in the field of one or more fixed monopoles. In the 
classical case, we found that this can be interpreted as geodesic motion on a manifold whose metric is the 
hyperbolic generalisation of the multi-centre metrics. We found that the Lagrangian for geodesic motion 




dx % dx'-' 



dr dr 



+ 1 -(e-V)E*= e -(E 2 -^). 



(4.49) 
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in the one-centre case admits an extra conserved vector of Runge-Lenz type, and we found the most 
general potential which could be added such that this persists. We shall solve the Schrodinger equation 
for motion in the one-centre metric with the potential found above. We shall find the bound state energies 
and also the scattering amplitudes for a particle scattered by the fixed centre. As special cases, we have 
the case of two monopoles, and also the MIC-Kepler and Kepler problems in hyperbolic space. We shall 
also show that the Schrodinger equation for the two centre case separates for pure geodesic motion and 
with a potential generalising that of the one-centre case. 

5.1 Bound states for the one-centre problem 

We wish to consider the quantum mechanics of a particle moving in a space whose metric is given by: 

ds 2 = Vh + V~ l (dr + uf . (5.1) 
With h the metric on H 3 . The Laplace operator of this metric may be written in the form ll42l : 

^ = + <«> 

where A h is the 'twisted' Laplacian on H 3 given by 



(1-r 2 ) 2 ( d 2 d\ l-r 2 ~ 
= —T^ [q-/ \ r ) + (5-3) 

where for the second equality we take Beltrami coordinates on H 3 and use the fact that uj r = to split 
the twisted Laplacian on H 3 into a radial part and a twisted Laplacian on S 2 . The time independent 
Schrodinger equation for motion on the manifold with metric (15.11) in a potential W is: 

'^V 2 + W 



Using (15.21 ) this becomes: 



A h + V 2 ^r - 2WV + 2EV )V = 0. (5.5) 







2 



dr 2 J 

We wish to consider a Hamiltonian which classically admits a Runge-Lenz vector, so we take 



V(r) = t + H, W(r) = ^ + ^ + ^(k + ^y 1 -^(k + ^^ (5.6) 

We will also impose the relation between V and u given by (12.231) . which implies that a = p. In order to 
separate variables for the Schrodinger equation, we set t = pip and we make the following ansatz for a 
solution: 

9 J KM (r,9,M) = R K M{r)D J KM {9A^), (5-7) 
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where there is no sum implied over repeated indices. The functions D J KM are the Wigner functions. We 
replace the angular derivatives by angular momentum operators according to 



^ = -(L\ + Ll), 

Acting on the Wigner functions we have that 

Ll + L 2 2 + Ll = J(J+1) 



-T 2 

-iv 3 . 



L\ = K 2 = eV- 



(5.8) 



(5.9) 



Since K is an integer or half-integer, we have a generalisation of the well known result of Dirac that 
the existence of a magnetic monopole implies charge quantisation once quantum mechanics is taken into 
account. We also have the usual relation that —J<K< J. The Schrodinger equation now reduces to 
an ordinary differential equation for the radial function: 



1 



r 2 ) 2 d 



d 



A B 



dr V dr 



where we have introduced some new constants 



R 



0), 



(5.10) 



A = J(J + 1) 
B = 2Ep-2(3 

C = 2Ek-2j-p 2 e 2 + J(J+1). (5.11) 
Equation (15.101 ) can be solved in terms of a hypergeometric function. The solution regular at r = is: 

2r \ J /i _ ~\ £(i+Vi+a-b-c) 



R J KM (r) 



where 



1 _ r \ 3^TVitfl-JJ-oj 2 r 

— — F(a,b;2J + 2;-—), 

1 + r J 1 + r 



^(1 + Vl+4:A + Vl + A- B- C-Vl + A + B- C 
-(1 + Vl + 44 + Vl + ^4 - B - C + Vl + A + 5 - C. 



(5.12) 



(5.13) 



In order that the solution is also regular at r = 1, we require that a is a non-positive integer. It is 
convenient to write a = —N + J + 1. This condition allows us to solve for E. Doing so, and replacing 
the radius of the hyperbolic space Rq by dimensional analysis we find that 



E 



N 



p(3 - kN 2 ± N^k 2 N 2 - 2kp(3 + - N 2 + e 2 p 2 ) + 2-?p 2 

p2 



(5.14) 



is a bound state energy level, provided that N takes one of the values \ep\ + 1, \ep\ + 2, . . . and also that 



-13 + pE N -N 2 >0 



(5.15) 
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is satisfied. We note a larger degeneracy than one might expect for the spectrum of this problem. The 
energy depends only on the quantum number N, not on the total angular momentum. This degeneracy is 
exactly analogous to the degeneracy in the spectrum of the hydrogen atom and is also due to the hidden 
symmetry whose classical manifestation is the Runge-Lenz vector. 

We may consider some of the special cases mentioned earlier. In the limit of pure geodesic motion, 
which describes the interaction of a hyperbolic monopole with a fixed monopole at the origin, we find 
that the bound states have energy given by: 




E N = ^ ( eV - N 2 ± N^N 2 - eV + - ^ 2 + eV) ) (5-16) 

with k — 1 + p/Ro- Taking the limit R — > oo, we recover the results for the energy levels of bound 
states of two Euclidean monopoles as derived in [16]. As in the Euclidean case, we discard the states cor- 
responding to the negative square root in (15.161) as these are tightly bound states which are not described 
by our well separated monopole assumptions. 

We can also recover the MICZ-Kepler problem on H 3 as a limit of of our system. This problem was 
studied by Kurochkin and Otchik [12J and by Meng in other dimensions ll43ll . The limit we require is 
the limit p — > 0, while keeping ep = p fixed. This only makes sense in our expression for En (15.141) if 
we take the positive sign for the square root. Taking this limit, we find the spectrum for the MIC-Kepler 
problem in hyperbolic space after setting 7 = — /i 2 /2 to remove an additive constant from the energy: 

$2 N 2 — I 

En = ~2N~ 2 ~ ~2RT' # = H + 1, H + 2, • • • , W-PRo\ , (5.17) 

where [x\ is the largest integer not greater than x. We may take the limit R — >• 00 and we get the 
spectrum for the MIC-Kepler problem in Euclidean space. Finally if we take Rq — > 00 and p = we 
have the spectrum of the hydrogen atom, for some appropriate choice of units. 

One interesting result of the above calculations is that the hydrogen atom in hyperbolic space only 
has a finite number of bound state energy levels, the number of levels being bounded above by ^—(3R . 
A curious fact of the hydrogen atom in flat space is that the partition function for a canonical ensemble 
of electrons at temperature T populating hydrogen atom energy levels does not converge: 



Z=Y^e^l TN \ (5.18) 



JV=1 



If the hydrogen atom is in hyperbolic space however, the partition function converges since the sum is 
cut off at some finite value of N. This is reminiscent of large black holes in AdS, which may be thought 
of as having been put in a 'box' in order to stabilise them against losing their energy through Hawking 
radiation. 



5.2 Schrodinger's equation in pseudoparabolic coordinates 

When considering Rutherford scattering in Euclidean space, i.e. the problem of an alpha particle scattered 
by a much heavier nucleus, it is convenient to work with parabolic coordinates, £,77. These are defined in 
terms of the usual radial coordinate r and the z coordinate by 

£ = r + z, rj = r — z. (5.19) 
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In these coordinates the Schrodinger equation separates and it is possible to solve the resulting ODEs to 
find the Rutherford scattering formula ll44ll . 

In hyperbolic space, the appropriate choice of coordinates for this scattering problem are the pseudo- 
parabolic coordinates defined in section lA.1.51 As in the previous section, we have that Schrodinger' s 
equation takes the form: 

Ah + y2 §^ -2WV + 2EV^j * = (5.20) 
with 

A h = h« (v, - uH^j (Vi - . (5.21) 

Now the coordinates on hyperbolic space are /i, v as defined in section lA.1.51 A very similar calculation 
to that performed in section EOI to separate variables for the Hamilton- Jacobi equation leads to separation 
of variables for the Schrodinger equation. More precisely, we can write 

(//, u, <f), t) = e im(t> e ieT M{fi)N(ij) (5.22) 

and the Schrodinger equation (15.201) implies the two second order ordinary differential equations 

1 + pi 2 d f^dM\ (m + ep) 2 M 2E(k-p) + 2f3 - 2 7 - e 2 p 2 ^ = RM 



pi dpi \ dpi J pi 2 1 + /i 2 

Wrf fJN\ _ (m -e P f N + 2E { k + P ) -20-2,- = _ KM (52J) 

v dv \ dv J v z 1 — v l 

K is here a separation constant. 

These equations may be put into the form of the hypergeometric equation and solved. The solution 
for M(pi), regular at pi = 0, is given in terms of the hypergeometric function by 

M(pi) = /il m+ep l(l + p?)-W m +zpHSK F ( at, h; Cl ; t^—] ■ (5.24) 
New constants have been defined: 

1 , 1/7-11 , T 

ai = - |m + ep\ + -VK +- + -V1-A 

1 1 r— 1 1 , T 

bt = -\m + ep\ + -VK -\ V 1 - A 

2 1 n 2 2 2 

Ci = 1 + \m + ep\ 

A = 2E(k - p) + 2f3 - 2 7 - e 2 /. (5.25) 
Similarly, the solution for N(u) regular at v = is given by 



N(u) = u lm ~ epl (l - u 2 )2+2^^F (a 2 ,b 2 ;c 2 ;v 2 ) (5.26) 

with 



1 11 / VK 

a 2 = -\m-ep\ + - + -VI- B - — 



1 11 , VK 

b 2 = -\m-ep\ + - + -Vl-B + — 

c 2 = 1 + \m — ep\ 

B = 2E(k + p) - 2p - 2 7 - e 2 p 2 . (5.27) 
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We can recover the bound state energy levels from these formulae by imposing that the wavefunction is 
also regular at fi — oo and v = 1. This will occur if: 

a\ = —rii, a 2 = — ri2 (5.28) 

where n 1 ,n 2 are integers called the parabolic quantum numbers. This gives exactly the same energies 
for the bound states as in the radial calculation and a counting of states shows that they have the same 
degeneracy. 



5.3 Scattering States 

We are more interested in the scattering states. In flat space, we seek to write the wavefunction at large r 
as the sum of an incoming plane wave and an outgoing scattered spherical wave: 



^ ^ e ~ ifcz + /(#)— , (5.29) 



we then interpret f{6) as the amplitude per unit solid angle for a particle to be scattered at angle 9 to the 
incoming wave. In hyperbolic space, the situation is very similar. Instead of plane waves, the appropriate 
objects in H 3 are horospherical waves which are constant on the horospheres (see section lAT2l . To discuss 
the scattering solutions of Schrodinger's equation, it is convenient to use the coordinates (x, 0, <p) for H 3 , 
in which the metric takes the form: 

h = d X 2 + sinh 2 X {dd 2 + sin 2 Odtf) . (5.30) 

In these coordinates, the horospherical wave solution of Schrodinger's equation in hyperbolic space with 
no potential is given by: 

\P = (cosh x - sinh x cos Qy l ~ iK ' (5.31) 

and the spherical wave solution is: 



* = — r— , (5.32) 
sinn X 

where k' is related to the energy E by k' = \J2E — 1. Thus we shall seek a solution to the full 
Schrodinger equation (15.201) which has the asymptotic form as x —> 00 

$ ~ (cosh x ~ sinh X cos Oy 1 -^' + f(ff)-^- (5.33) 

smb. x 

and we will again interpret f{6) as a scattering amplitude. In order to find a solution with these asymp- 
totics, we look once again at the solution in pseudoparabolic coordinates. It is straightforward to show 
that (15.311) in pseudoparabolic coordinates becomes 

1 IK,' 

[(l+/i 2 )(l-i/ 2 )p + ^. (5.34) 

The condition on the asymptotic form of the wavefunction as x ^ oo becomes a condition as v — > 1 . 
We define k and A to be 



k = V 2E (P + k)-l-2j-2p + e y, A = v / 2E(k - p) - 1 - 2 7 + 2fJ - e 2 p 2 . (5.35) 
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In order to get the correct asymptotics, we set m = —ep and K — — 1 — iX. The solutions (15 .24115^61) 
then simplify and we get: 

M{n) = (l + /i 2 )^ 

N(y) = ^ 2|ep| (l-z/ 2 )HfF^|ep| + l + ^(fi: + A),|ep| + ^( ft :-A);l + 2|ep|;z/ 2 ^ (5.36) 
As v — > 1, we can expand the hypergeometic function using the formula: 

F(a , b; c; X) *Z' ^-"-g + (1 - «)— r(c) ^r C) . (»7, 

r(c — a)r(c — o) r(a)r(o) 

Performing the expansion and multiplying by a factor to normalise the incoming wave, we find 

"~ 1 [(1 + /i 2 )(l - ^ 2 )] (1 + /i 2 )^ A " K ) + A(l - z/ 2 )5"f (1 + /i 2 )^, (5.38) 

where 

VttxW (\ P -n\ - Uk 4- \\\ r fUnl 4- 1 - if* - XYl 

(5.39) 



r(m)r ( ep 


-i(« + A))r(|ep 


+ 1 • 


-I(k-a)) 


r (—in) r ( ej 


»| + 1 + |(« + A))r(|ep 


+K«-a)) 



We can now return to the (%, 0, 0) coordinates, and we find that this formula becomes: 

* x ~°° (cosh X ~ sinh X cos 0)-i-« e -<(A-«) lnsin f + _ CO sec 2 - e - i(A - K) lnsin * — — (5.40) 

2 2 sinh % 

We thus find the scattering amplitude to be: 

f{6) = ycosec 2 ^. (5.41) 

It is interesting to note that both the incoming and outgoing wave are distorted relative to the proposed 
asymptotic form (I5.33I) . This occurs also in the Euclidean case, because the Coulomb potential does not 
decay fast enough to justify an expansion in terms of free waves far from the scattering centre. 

It is possible to recover the bound state energy levels by making the analytic continuations k — > ik 
and A — > i\. We then find that the scattering amplitude has poles at precisely the values of E which 
were found previously to be bound state energies. Another useful check of this formula is to use it to 
find the Rutherford scattering formula. We first set p = 7 = e = and k = j3 = 1, then we replace the 
hyperbolic radius R by dimensional analysis. Finally we let R — > 00 and we find 

f(ffs _ 1 r (! + W) r c 4r> 

I[ ) 2fc' 2 sin 2 f r(l-V^) } 

where E = \k' 2 . This is precisely Rutherford's formula, as found in ll4*4l . 
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6 The two centre problem 



6.1 Separability of Schrodinger 's equation for the two-centre problem 

In P2l it was shown that the Schrodinger equation for the two-centre Taub-NUT metric separates in 
spheroidal coordinates. This generalises the fact that the Schrodinger equation for a diatomic molecule 
separates in the same way that the Runge-Lenz vector for Taub-NUT generalises that of the Hydrogen 
atom. We find here that in a suitable set of coordinates, the two-centre LeBrun metric has a separable 
Schrodinger equation. The coordinates are the pseudo spheroidal coordinates used by Vozmischeva P31. 
which for convenience are also described in section (IA.1.41) 
We consider motion in a space with metric 

ds 2 = Vh + V~ 1 {dT + uf (6.1) 

Where V and to are given by 

V(P) = k- pi cothP (P, Pi) - P2 cothD (P,P 2 ), du = * h dV. (6.2) 

It is convenient to work for the moment with the pseudosphere model of Hyperbolic space (see section 
IA. 1.11) . The general point P G H 3 has coordinates (W, X, Y, Z) G E 3 ' 1 which are subject to the constraint 
W 2 — X 2 — Y 2 — Z 2 = 1. We suppose that the fixed centres are at (a, 0, 0, ±/5). The distance between 
two points P, P' is given in terms of the 5*0(3, 1) invariant inner product on E 3,1 by: 

(P P') 

D(P, P') = 9, where cosh 9 = n^jlT^lj = WW> ~ XX ' ~ YY ' ~ ZZ '- ( 6 - 3 ) 
Thus the potential function V takes the form in these coordinates 

V = k - pi coth 8+ - p 2 coth 0_ , with cosh 6± = aW =F (3Z. (6.4) 

In order to pass to the pseudospheroidal coordinates, we first take polar coordinates (P, <p) in the X-Y 
plane. The pseudospheroidal coordinates £, 77 are then defined in terms of W, P, Z by the relations (IA-121) : 



W = -^(a 2 -e)(a 2 + V 
a 

Z = S -^jp-y/(P 2 -e)(P 2 + V 2 ), 



,2> 



p = §v_ 

The metric on H 3 in these coordinates is given by (IA- 1 81) 



(6.5) 



Using these relations with (16.41) it is possible after some algebra to find the following formula for V in 
terms of the pseudospheroidal coordinates: 



V = k- ( Pl + P2 ) sgn(Z)(p! - p 2 ) -^-^ . 
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The simplest way to compute cu is directly from the formula in ( 16.21) . Using ( I6.6I ). again after some 
algebra, it can be shown that u takes the form: 

^ = ^3 I (Pl +:P2)r/ sgn(Z)( Pl -p 2 )e I #. (6.8) 

We consider the Schrodinger equation for the metric d6.ll) with a potential W. The Schrodinger 
equation can be re- written in the form (15. 201) . We make a separation ansatz 

^ me (£, V, r, 0) = e im ^ eT X(0Y(v) (6-9) 
and we find that the Schrodinger equation separates when W satisfies 

W = 7 + /3+coth0 + + /3_coth0_. (6.10) 
In this case, the functions X(£) and Y(rj) obey the ordinary differential equations 

KX = - y - a | _ ew _ a «) + ^ _ ^ x 

+ (Vsgn(Z) + ^Pj v/( a 2-e)(/5 2 -e)^ (6-11) 

and 

-KY = +riW + ri) ^) + \ V ~^f) 

+ (c~ + sgn(Z)^ y/(a 2 + rj 2 )(0 2 + ri 2 )Y. (6.12) 

is a separation constant. In the case where the potential vanishes, K corresponds to a Killing tensor of 
the metric (16.11) . We have defined some new constants 

A ± = 2Ek-2-f - e 2 k 2 - e 2 (p 1 TP2) 2 , 

B± = a 2 (3 2 (m 2 + e 2 (p 1 ±p 2 ) 2 ) J 

C± = 2e 2 ( Pl TP2)k-2E( Pl ±p 2 )-2(f3 1 TP2), 

D ± = ±2emap(p 1 ±p 2 ). (6.13) 



6.2 Separability of the Hamilton- Jacobi equation for the two-centre problem 

Additive separability of the Hamilton- Jacobi equation follows, formally at least, from the multiplicative 
separability of the Schrodinger equation using a semi-classical approximation. It is however straight- 
forward, using the expressions for V and u to separate variables for the Hamilton- Jacobi equation in 
pseudo spheroidal coordinates explicitly. We find that the Hamilton- Jacobi function takes the form: 

S = ^ + rp T + M(£) + N(rj), (6.14) 
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where M, N satisfy ordinary differential equations: 



dM\ 2 ^ b + 



+ (Vsgn(Z) + ^ ^/( a * (6-15) 



and 



K = (« 2 + ^)(^1 2 )(f) 2 + 2 + ^ 

+ (c~ + sgn(Z)^ y/^ + rj^iP + rj 2 ). (6.16) 



The new constants are 



«± = p T 2 +p T 2 (p lT p 2 ) 2 + 2 7 -2£/ ; ; 2 , 

b ± = a 2 (3 2 (p/ + ( Pl ±p 2 ) 2 p T 2 ), 

c* = 2£(pi=Fp 2 ) +2(^=F/3 2 ) -2kp T 2 ( Pl TP2), 

d ± = ±2af3p <j) p T {p l ±p 2 ). (6.17) 

Thus the two centre problem, with a potential generalised from that of the one centre problem, is 
classically integrable in the Liouville sense, i.e. there are sufficiently many conserved quantities to reduce 
the problem to quadratures. The 'extra' conserved quantity, K corresponds in the absence of a potential 
to a rank 2 Killing tensor of the metric (16.11) . 

The system is also quantum mechanically integrable. In the operator formalism, this corresponds to 
being able to find a basis of states which are eigenstates of operators which commute with the quantum 
Hamiltonian. This quantum integrability is a stronger result than classical integrability. In ll46l Carter 
showed that a classical rank 2 Killing tensor will not in general generate a conserved charge in the 
quantum regime. There is an anomaly given by: 



K,H 



oc 



For a Ricci flat background, this anomaly will always vanish and there will be a conserved quantum 
charge associated to the classical conserved quantity. In our case the Ricci tensor does not vanish so a 
priori we cannot expect a classical symmetry to give rise to a quantum conserved charge. The separation 
of variables for the Schrodinger equation shows however that the anomaly vanishes in this case. 



7 Conclusion 

We have investigated the motion of hyperbolic monopoles in the limit where the monopoles are well 
separated and found many similarities with monopoles in Euclidean space. Because of the absence of 
boost symmetries for hyperbolic space it was necessary to consider the case where N monopoles are 
fixed while one is free to move in the asymptotic field of the others. We found that this motion may 
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be interpreted as a geodesic motion on a space whose metric is a LeBrun metric of negative mass. We 
interpret this as the metric on a non-geodesic surface contained within the full moduli space of the theory. 

We then investigated in depth the classical and quantum motion of a monopole in the presence of a 
single fixed monopole at the origin. It was found that this system may be considered within a broader 
class of systems which have similar properties arising from the existence of an extra 'hidden' constant of 
the motion which generalises the Runge-Lenz vector of the classical Kepler problem. These systems have 
an enhanced dynamical symmetry algebra and this results in some attractive properties. The system is 
integrable in the Liouville sense and one can show that the orbits, when expressed in suitable coordinates, 
are conic sections. The hodograph, suitably defined, is elliptic and becomes circular precisely when the 
system under consideration is the hyperbolic Kepler system. The Hamilton-Jacobi equation separates 
in both spherically symmetric coordinates and pseudoparabolic coordinates, enabling the surfaces of 
constant action for scattering particles to be found analytically. 

The classical integrability persists in the quantum regime and the Schrodinger equation separates in 
spherically symmetric and pseudoparabolic coordinates. This enables us to calculate the bound state 
energies and the scattering amplitude for a monopole moving in the field of another fixed at the origin. 
We find similarities with the hydrogen atom problem: the degeneracy of the energy levels is enhanced by 
the existence of the Runge-Lenz vector; the scattering amplitudes have the same angular dependence as 
appears in the Rutherford formula. 

We have also considered the problem of a monopole moving in the field of two fixed monopoles and 
we find that it is integrable both as a classical and a quantum system. Geometrically the integrability 
properties of both the one and two centre cases are due to the existence of rank 2 Killing tensors. These 
tensors generate non-geometric transformations of the phase space which are symmetries of the system 
and give rise to conserved quantities which are quadratic in the momenta. 

Thus we have found many similarities between the motion of a hyperbolic monopole about a fixed 
monopole and the same problem in flat space. The systems we considered generalise the original Kepler 
problem and exhibit many of beautiful properties which arise as a consequence of its symmetries. 
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A Some useful properties of Hyperbolic Space 



A.l Coordinates on Hyperbolic space 

For convenience, we collect here several coordinate systems for hyperbolic space which will prove useful 
in this work. 

A.l.l The pseudosphere model 

The most easily visualised model of hyperbolic space is as the upper leaf of the unit pseudosphere in 
Minkowski space. This model shows most clearly the parallels between spherical and hyperbolic geom- 
etry: 

H 3 = {(W,X,Y,Z) E E 3,1 : W 2 -X 2 - Y 2 - Z 2 = 1,W > 0} (A-l) 
where the metric on E 3,1 is given by 

ds 2 = -dW 2 + dX 2 + dY 2 + dZ 2 . (A-2) 

The metric on H 3 is then the restriction of this metric to the tangent space of the pseudosphere. Lorentz 
transformations of the full space E 3,1 preserve both the metric (IA-21) and the pseudosphere condition 
(IA-11) and so correspond to isometries of H 3 . We see that SO (3, 1) acts transitively on H 3 , and we may 
write H 3 = SO(3, l)/50(3). As in the case of the three sphere in E 4 , the geodesies are given by the 
intersections of 2-planes through the origin with the pseudosphere. It is useful to define the various 
coordinate systems that we will use in terms of this model. 

A.1.2 Poincare Coordinates 

Among the best known of the coordinate systems on H 3 are the Poincare coordinates. These are obtained 
from the pseudosphere model by the analogue of stereographic projection. 



H.,_----r4(W, 0,0,0) 





(1, 0, 0, 0) 








'(-1,0,0,0) 



Figure 3: Stereographic projection of H 3 onto plane through the origin 
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We consider the projection through the point (—1, 0, 0, 0) onto the plane with normal (1, 0, 0, 0) in 
E 3,1 which passes through the origin. This situation is shown with one dimension suppressed in Figure 
13 From the diagram, we see that the point P in the pseudosphere is mapped to a point P' on the plane 
W = 0. We introduce cylindrical polar coordinates for E 3,1 according to: 

X = R sin 6 sin $ 
Y = R sin cos $ 

Z = RcosS. (A-3) 

We denote the coordinates of P' by (0, r, 9, <p). Clearly from the diagram we have that the coordinates of 
P' are related to those of P by 



R r 



W + 1 1 



9 = 0, $ = <j). (A-4) 



Where the condition that P lie on the pseudosphere becomes W 2 — R 2 = 1. We note that since W > R, 
we must have r < 1, so H 3 has been mapped to the interior of the unit ball in M 3 . The sphere r = 1 
corresponds to the 'sphere at infinity' of H 3 . In the r, 9, <p coordinates the metric on H 3 can be written as 

ds 2 = _ 4 ^ 2 {dr 2 + r 2 (d9 2 + sin 2 9d0 2 )} . (A-5) 

This is the model of H 3 known as the Poincare ball. We recognise the term in braces as the flat metric on 
IR 3 so that the metric is conformally flat in these coordinates. The geodesies in this model correspond to 
the arcs of circles which intersect the sphere r = 1 orthogonally. 



A.1.3 Beltrami Coordinates 

The Beltrami coordinates on H 3 are the hyperbolic analogue of the gnomonic coordinates on S 3 . They 
are again defined by a projection, but this time we project through the origin onto the tangent plane to H 3 
at (1, 0, 0, 0), as shown in FigurelU 




Figure 4: Gnomonic projection of H 3 onto the tangent plane at (1, 0, 0, 0) 

We can again write the coordinates of P' in terms of those of P. We note that the projection will 
change the magnitude, but not the direction of the vector X = (X, Y, Z). This is simply the observation 
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that the angles $ and are unchanged by the projection. The Beltrami coordinates r are defined in terms 
of the coordinates in E 3 ' 1 according to: 

r = — (A-6) 

W 

Once again, the condition that P lies on the pseudosphere is W 2 — X 2 = 1, so we are again mapping H 3 
to the interior of the unit ball. The metric in Beltrami coordinates is given by 



_ dv 2 (r ■ drf 

1 i 1 2 / , ,9x2- V A ') 

l-|r| (l-|r|) 



Defining spherical polar coordinates on M 3 in the usual way, the metric may alternatively be written 

dr 2 r 2 

ds 2 = — 5 + (d9 2 + sin 2 6d<p 2 ) (A-8) 

(1 -r 2 ) 2 1 ~r 2( ' 

Since the geodesies on H 3 are given by the intersection with 2-planes through the origin in the pseudo- 
sphere model, it is straightforward to see that they project back to straight lines in the Beltrami ball 
model. 



A.1.4 Pseudospheroidal Coordinates 

When considering the two centre problem in flat space, it is convenient to work in spheroidal coordinates. 
Since the two centres may be taken to lie at (0, 0, ±a), we take the polar angle as one of the coordinates. 
The other two coordinates are defined in terms of the distances to the centres, r + and r_ by 

C = ^±^, A=^±^. (A-9) 
s 2a 2a 

The coordinate surfaces £ = const, are prolate spheroids and the surfaces A = const, are hyperboloids. 
In these coordinates the Laplacian on flat space separates. 

In our investigations we shall require an analogous set of coordinates for H 3 . We shall use the coor- 
dinates defined by Vozmischeva in P31. We first set 

X = Pcos0, F = Psin0, (A- 10) 

so that the metric on E 3,1 becomes 

ds 2 = -dW 2 + dZ 2 + dP 2 + P 2 d<p 2 (A- 1 1 ) 

and the equation of the pseudosphere becomes W 2 — P 2 — Z 2 = 1. We assume that we have two centres 
located at (a, 0, 0, ±(3), where of course a 2 — f3 2 = 1. The pseudospheroidal coordinates £, 77 are defined 
by 



W 



-V(a 2 -e)(a 2 + V 2 ), 
a 



P = % (A-12) 
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with 

0<£</3, < r/ < oo (A-13) 

Since sgn(Z) enters these equations, the equations (IA-121) in fact define two coordinate patches, which 
together cover H 3 and overlap along the plane Z = 0. From these equations, we can deduce the relations 



W 7 



p2 

£2 : 



a 2 + 7] 2 P 2 + rf' 



+ 



if 



(A-14) 



In order to interpret the geometric meaning of the £, 77 coordinates, we may use the gnomonic projection 
and consider the coordinate surfaces in Beltrami coordinates. Using cylindrical polars on the Beltrami 
ball, from above we have that 

Z P 

Z = W> p= w> (A " 15) 

while (p remains unchanged. Thus, in Beltrami coordinates, the equations which determine the coordinate 
surfaces of £ and rj become: 



1 



ft 



ft" 



p 2 -e 

z 2 
(3 2 + ri< 



El 

e 



+ 



(A-16) 
(A-17) 




Figure 5: Coordinate curves £ = const, and r] = const. 



As £ and r\ vary, these describe a family of hyperbolae and ellipses respectively. A plot of the coordi- 
nate curves in the z — p plane is shown in Figure where for convenience we allow p to be negative. 

Using equations (IA-111) and (IA-121) . one can show that the metric on H 3 in pseudospheroidal coordi- 
nates is 

*'= v <A - 18) 
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A.1.5 Pseudoparabolic Coordinates 



For one-centre problems in flat space with a distinguished spatial direction, it is often useful to use 
parabolic coordinates. These may be thought of as spheroidal coordinates where one of the two centres 
has been allowed to recede to infinity, while the other is fixed at the origin. Again following Vozmischeva 
P31 . we shall construct an analogous set of coordinates on H 3 from the pseudospheroidal coordinates 
given above. 

We first apply a Lorentz transformation to E 3 ' 1 in order to move the centre at (a, 0, 0, —j3) to (1, 0, 0, 0). 
Without loss of generality we may write a = cosh & j3 — sinh & The point at (W, Z, P, 0) moves to 
(W, Z', P' , <//), where the new coordinates are related to the old by: 

W = cosh -W - sinh -Z', 
2 2 

Z = cosh ^Z'- sinh | W, 

P = P', = 0'. (A- 19) 

We also rescale the variables £ and rj by 

£ = a£' = cosh — rj = (3r]' = sinh —7/ (A-20) 

In order to find the new coordinate surfaces, we substitute (IA-191) and (I A- 201) into (IA-141) and take the 
limit as ip — > 00, keeping terms of order e~'^' . The new coordinate surfaces are then given by: 

2Z'(W'-Z') = J^(W'-Z') 2 - l -^^P' 2 = ^\W'-Z') 2 -^ (A-21) 

2Z'(W'-Z') = !^(W'-Z') 2 + ^^P' 2 = -v 2 (W'-Z') 2 + ^- (A-22) 

1 + 7] 2 Tj V 1 

Where we have defined new coordinates [i = and v = j^p- We now have no further use for 
the unprimed coordinates, so drop the primes. We can solve equations (IA-211) . (IA-221) . together with the 
pseudosphere constraint W 2 — P 2 — Z 2 = 1 to express the points on the pseudosphere in E 31 in terms 
of \i and v. This gives 

\y — 2 + ^~^ 2 p _ ^ z _ V 2 -" 2 (A _ 23) 

2^(1-^(1+^2)' y/(l-V*)(l+fy 2^(1-^(1+^2)' 

where 

0</i<cx), < v < 1. (A-24) 

This time the coordinates cover the whole of H 3 . The sphere at infinity is given by v = 1. We may 
once again use the gnomonic projection in order to visualise the coordinate surfaces as surfaces within 
the Beltrami ball model. We find that the Beltrami coordinates z, p are given in terms of /x, v by: 



Z n 2 - v 2 P _ pv 

W ~ 2 + /i 2 - z/ 2 ' 9 ~ W ~ 2 + /j, 2 - v 



P = 777 = o ■ " 2 _„ 2 - (A-25) 
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Figure 6: Coordinate curves fx = const, and v = const. 



Either from these equations, or from the projected versions of equations (IA-211) . (IA-221) . we find that the 
surfaces fj, = const, are hyperboloids and those of v = const, are ellipsoids. The plane 4> = is shown 
in Figure |6j 

Using equations (IA-231) and (IA-1 lb . it is a matter of straightforward calculation to find the metric on 
H 3 in the /i, v, coordinates. We find that 



da* 



/i 2 + V 2 



dfi 2 



dv 2 



(A-26) 



A.2 Horospheres 

When considering scattering processes in Euclidean space, the incoming wave is usually considered to be 
a plane wave, i.e. a wave whose phase is constant along any plane parallel to some vector. The important 
property of such planes is that they are everywhere normal to a set of parallel geodesies. Whether we 
are considering classical action waves or wavefunctions, the interpretation is that particles are moving 
parallel to some fixed direction. For hyperbolic space, we would like a similar set of waves to describe 
incoming particles initially parallel to some direction. In order to do this, we consider a pencil of parallel 
geodesies in H 3 , which meet at a point on dM 3 . The horospheres are a set of surfaces everywhere normal 
to these geodesies. It is most convenient to visualise these surfaces in the Poincare ball model, since this 
is conformally flat so angles are as one would expect in Euclidean space. A set of parallel geodesies in 
this model are a set of circular arcs which meet the sphere r = 1 orthogonally at some given point, say 
(0, 0, 1). A surface which is everywhere normal to this pencil of geodesies is a sphere which is tangent 
to the sphere r = 1 at the point (0, 0, 1). A picture of the y = plane is shown in figure|7] 

The horospheres may be thought of as a limit of a set of spheres in H 3 whose centre goes to infinity 
whilst the radius also tends to infinity. The horospheres carry a natural Euclidean geometry and the 
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Figure 7: A plot showing horospheres and their normal geodesies in Poincare coordinates 



restriction of the hyperbolic metric to each horosphere is flat, a result known as Wachter's theorem H71 
l48l . Clearly there is nothing special about our choice of boundary point, so there is a set of horospheres 
associated with every point on S 2 ^ = <9H 3 . 

A brief calculation shows that the horospheres associated with a pencil of parallel geodesies originat- 
ing from (0, 0, 1) are given in pseudoparabolic coordinates by the equation 

(1 + /i 2 )(l - v 2 ) = const. (A-27) 
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